Quantifying Entanglement in Quantum Complex Networks 
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We introduce the use of entanglement entropy as a tool for studying the amount of information 
stored in quantum complex networks. By considering the ground state of a network of coupled 
quantum harmonic oscillators, we compute the information that each node has on the rest of the 
system. We show both analytically and numerically that the nodes storing the largest amount of 
information are not the ones with the highest connectivity, but those with intermediate connectivity 
thus breaking down the usual hierarchical picture of classical networks. As a byproduct, our results 
point out that the amount of information available for an external controller connecting to a quantum 
network is bounded although it can be maximized by means of a moderate number of connections. 

PACS numbers: 89.75. Fb,03.67.-a,89.70.Cf 



The advent of network science has influenced the re- 
search in many fields of science ingeneral, and physics 
in particular, in a pervasive way Since the discov- 
ery of the structural features of real social, biological and 
technological networks @, [|[ , the development of the the- 
oretical machinery of network science has blossomed as 
an efficient framework to interpret the many interaction 
patterns encoded in real-scale complex systems of diverse 
nature [H and to model correctlythe dynamical processes 
taking place on top of them pll(| . 

One of the most important avenues of research in net- 
work science is its connection with information theory. 
In this way, different information-theoretical tools have 
been proposed to characterize the complexity of networks 
beyond the typical statistical indicators such as their de- 
gree distribution, clustering coefficient, degree correla- 
tions, etc For instance, Shannon entropy, as shown in 
[7H10l|. has been successfully applied to characterize the 
complexity of ensembles of networks sharing some struc- 
tural features while information-theoretical tools have 
been also applied to the study of diffusion processes on 
top of networks, such as random walks [lll - [l6j . 

The synergy between the field of complex networks and 
that of information theory has recently appealed to the 
quantum information community. As a p roduct, clas- 
sical results on percolation theory |17h2(J and network 
science, such as the small- world effect [2 ij . have been re- 
visited in networked structures of coupled quantum sys- 
tems as a first step for designing quantum communica- 
tion networks. Conversely, the use of quantum dynamical 
processes, such as quantum random walks |22| and their 
application to rank the importance of network elements 
[23M26I I , has given new quantum information perspectives 
to classical problems of the network realm. 

In this Letter we quantify the amount of information 



that a single element of a quantum network shares with 
the rest of the system. To this aim we consider a network 
of quantum harmonic oscillators and analyze its ground 
state to compute the entropy of entanglement that vac- 
uum fluctuations creates between single nodes and the 
rest of the network. By using the Von Neumann entropy 
27, 28] we show that the amount of information available 
to nodes is limited, regardless of how large the number 
of connections of a node is. Moreover, in contrast to the 
usual scenario in classical networks, the nodes with the 
largest entanglement are not the hubs, but those com- 
ponents with intermediate connectivity. As a byproduct, 
our results show that in order to extract the maximum 
information of a quantum network it is enough to add 
an external controller connected to a moderate number 
of elements. Alternatively, a perturbation via quantum 
operations of intermediate-degree nodes induces the max- 
imum damage to the quantum network state. Finally we 
round off the Letter by commenting on further implica- 
tions of our results in quantum gravity models. 

Theory and Model.- As usual, we define a network as a 
set of N nodes and E edges (or links) accounting for their 
pairwise interactions. The network backbone is usually 
encoded in the Adjacency matrix, A, such that Ay = 1 
if an edge connects nodes i and j while Ay = other- 
wise. In this Letter we restrict to undirected networks 
so that Ay — Aji. Although matrix A stores all the 
structural meaning of a network it is more convenient to 
rely on the so-called network Laplacian, L, to analyze its 
structural and dynamical properties [231. The Laplacian 
of a network is defined from the Adjacency matrix as 
Lij = kiSij — Aij, where fc$ = Ay is the connectivity 
of node i, i.e., the number of nodes connected to i. 

Once defined the physical skeleton of the network we 
define its dynamical setup. As anticipated above, we con- 
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sider nodes as identical quantum oscillators, interacting 
as dictated by the network topology encoded in L. The 
Hamiltonian of the quantum network thus reads: 

^network = r (P T P + X T (I + 2 C L)x) , (1) 

here I is the N x N identity matrix, c is the cou- 
pling strength between connected oscillators while p T = 
(Pi,P2,—,Pn) and x T = (xx, ,x%, ...,xn) are the oper- 
ators corresponding to the momenta and positions of 
nodes respectively, satisfying the usual commutation re- 
lations: [x, p T ] — ihl (we set ft = 1 in the following). 
The characterization of the quantum dynamical states is 
given by the eigenvectors of the Laplacian matrix and 
their corresponding eigenvalues shifted by 1. Note that, 
although the classical ground state of the above system is 
trivial (all the particles having xi = Pi = 0), the quantum 
nature of the system prevents nodes from being at rest 
due to Heisenberg uncertainty principle. Moreover, the 
nature of the associated quantum vacuum fluctuations 
only depends on the pattern of physical interactions, i.e., 
the network structure. 

We now study the properties of the ground state to 
quantify the amount of information each of the elements 
of a network shares with the rest of the system via quan- 
tum fluctuations. To this aim, we consider the partition 
of the network into a node, say i, and its complement 
i c , i.e. the rest of the network. Then, we compute the 
mutual information shared by the two parties as: 

l(i\i c ) = Si + S lC - S tot . (2) 

Here Si and S& are marginal entropies and Stot is the 
total entropy of the network. It is natural to choose the 
Von Neumann entropy to quantify the quantum informa- 
tion of the system, yielding Stot = for the ground state 
(as it is a pure state). The marginal entropies for i and 
i c read [I|: 

Si = = (p-i+jj lo s(pi+^) -fa- o) lo g(i"i-^) > ( 3 ) 

which is a monotonically increasing function of /Ltj that is 
characterized by the second moments of the positions and 
momenta of nodes, /i, = \J (x\)(p 2 i ). After some algebra 
(see Appendix [A| we are able to quantify the value of /Ji 
as: 
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where {Xj} are the eigenvalues of the network Lapla- 
cian L and matrix S accounts for the normal mode 
transformation that diagonalizes the network Laplacian: 
L d = S T LS with S T S = I. 

From Eq.Q it is clear that each node has some mutual 
information with the rest of the system provided fii > 1/2 
whereas from Eq.Q we conclude that the amount of 
information depends on its contribution to each of the 
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FIG. 1. (color online) Average entropy (S) versus the coupling 
strength, c, for Random Regular Graphs (RRG) •, Erdos- 
Renyi (ER) ■ networks, and Configurational (SF-CONF) ♦ 
and Barabasi-Albert Scale-Free (SF-BA) A networks. The 
average connectivity of the networks is (k) — 4. 



Laplacian eigenvectors. We finish by noting that, since 
the total network is in its ground (and pure) state we have 
Si = Sic = I(i\i c )/2. Therefore, the information that a 
node shares with the network is intrinsically due to quan- 
tum correlations. Equivalently, the mutual information 
is, itself, a measure of the entanglement (quantified by 
Si) between a single node and the rest of the system. 

In the following we quantify the entanglement en- 
tropies of nodes embedded in different network topolo- 
gies. First, we explore two homogeneous network sub- 
strates: (i) Random Regular Graphs (RRG), in which 
all the nodes have the same number of contac ts ( kj = 
(fe),Vi), and (ii) Erdos-Renyi (ER) networks [31(, for 
which the probability of finding a node with k neighbors, 
P(k), follows a Poisson distribution so that most of the 
nodes have a degree k close to the average (k). Besides, 
we have analyzed networks having a scale-free (SF) pat- 
tern for the probability distribution, P(k) ~ fc~ 3 , con- 
structed by means of a configurational random model 
g-CONF) [32| and the Barabasi-Albert model (SF-BA) 

Macroscopic analysis.- First we report the average en- 
tanglement entropy per node (Si) — Yl%=i Si/N, as a 
function of the coupling between oscillators. From Fig. Q] 
it is clear that entanglement increases with the coupling 
strength c. However, fixing the value of the coupling, 
the entropy per node of each network shows differences. 
In particular, as the heterogeneity of the degree distri- 
bution P(k) of the network increases (from RRG to SF 
networks) the value of (Si) decreases. Thus, the presence 
of different connectivity classes within the same network 
seems to be detrimental for the (average) entanglement 
within its components. 

Microscopic analysis. - The above finding motivates the 
study of the entanglement entropy at the level of single 
nodes in networks having different connectivity classes 
such as ER and SF networks. To this aim, we collect 
the entanglement entropies of the Nk nodes having con- 
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FIG. 2. (color online) In (a) we show the microscopic picture: 
a particular node and its boundary (first neighbors) with the 
rest of the system (blurred). The rest of the panels show the 
average entropy of nodes with degree k, (Sk) for the following 
topologies: (b) ER, (c) SF-CONF and (d) SF-BA networks. 
Each symbol refers to a value of the coupling strength c. In 
particular, we have that • ■ ♦ A correspond to c = 0.2, 0.4, 
0.6 and 0.8 respectively. All the networks have the same av- 
erage degree (k) = 4. Solid lines represent the theoretical 
curves calculated using the mean-field formulation. 



nectivity k and define the average entanglement of the 
degree class k as: (Sk) = J2i\k i= k Si/Nk- One would ex- 
pect that the larger the connectivity k of a node the more 
correlated with the rest of the network it is, and thus the 
larger the value of (Sk) is. 

The growth of (Sk) with k would point out the exis- 
tence of an area law [3J] for complex networks. The area 
law states that the entanglement entropy of an inner re- 
gion of a system scales with the size of the boundary 
connecting it with the rest of the system. For a regular 
lattice, the size of the boundary of an element is given by 
twice its dimensionality thus, in analogy, for a node in a 
complex network its boundary [see Fig. [2] (a)] is given by 
its connectivity, i.e., (Sk) ~ k. Taking this picture, one 
can consider a heterogeneous network as a collection of 
different dimensionalities within the same system. 

The panels in Figs. [5] and [3] summarize our findings 
for the behavior of (Sk)- In Fig, [5] we explore ER [panel 
(b)], SF-CONF [panel (c)] and SF-BA [panel (d)] net- 
works for different values of the coupling c = 0.2, 0.4, 0.6 
and 0.8. For ER networks we observe that the value of 
(Sk) increases with k although we note that the growing 
trend seems to saturate for large values of k pointing out 
that entanglement is bounded. On the other hand for SF 
networks, displaying a larger heterogeneity for the collec- 
tion of degrees, the growing trend of (Sk) only holds for 
small to moderate connectivities fc, then (Sk) reaches a 
maximum and starts to decrease. Eventually, those nodes 
with sufficiently large k would drop its entanglement en- 
tropy. As a result, hubs are not the most entangled nodes, 
but there is an optimally-correlated class of nodes having 




FIG. 3. (color online) Average entropy of nodes with degree 
k, (Sk), for all the topologies under study: (a) RRG , (b) ER, 
(c) SF-CONF and (d) SF-BA networks. All the oscillators 
are coupled with the same strength c = 0.6. Each symbol 
represents a different value of the average degree (k) of the 
system. Solid lines represent the theoretical behavior calcu- 
lated with the mean-field formulation. 



moderate connectivity. The plots in Fig. [3] confirm the 
above results. In these cases, we have fixed c = 0.6 and 
changed the mean connectivity of the RRG [panel (a)], 
ER [panel (b)], SF-CONF [panel (c)] and SF-BA [panel 
(d)] networks |35j . 

Mean-field formulation.- In the following we develop a 
minimal model aimed at capturing the rise-and-fall be- 
havior of (Sk) in SF networks. We propose a generalized 
mean field approximation that goes beyond the standard 
one in which the many-body Hamiltonian is reduced by 
replacing the interaction terms by their mean value, thus 
hindering the effect of the heterogeneity in the number 
of inputs (degree) that nodes receive. 

The simplest framework for dealing with bipartite en- 
tanglement is shown in Fig. [2j (a) . We consider both the 
node of interest (say 0) and its first k neighbors, whereas 
the influence of the rest of the network (blurred in the 
figure) is taken into account by renormalizing the fre- 
quency of the latter ones. In this approximation, the 
local Hamiltonian of a node of degree k reads, 



%f — 2 ) Po 



where the renormalized frequency of the k neighbbors 
reads: v~ = 1 + c 2 k, where n is a fitting parameter (see 
below) . 

This model is exactly solvable (see Appendix |B|) . 
Therefore, we can find the entropy of the central node 
(Sk) MF analytically as a function of its connectivity k 
and k. In Figs. [2] and |3J we plot with solid lines the 
curves (Sk) MF obtained after tuning the single parame- 
ter k for each of the curves. Interestingly, the values of k 
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FIG. 4. (color online) Controlled entanglement via an exter- 
nal (blue) node coupled to an ER graph. In the plot we show 
(filled dots) the entanglement (Sk) of the nodes of the ER 
network as a function of their connectivity. In addition, we 
show (empty squares) the evolution of the entanglement of 
the external controller as a function of its connectivity, i.e. 
the number of links launched to the target network. 

used for each curve only depends on the network topol- 
ogy. As reported in Appendix [C] its value is independent 
of the value of c but proportional to the average conec- 
tivity (k). From the figures it is clear that our mean-field 
approximation agrees fairly good with numerics. 

In addition to the quantitative agreement, the analyt- 
ical estimation (Sk) MF allows to explain the rise-and- 
fall of entropy across degree classes. As shown in the 
Appendix [C] this phenomenon lies in the fact that hubs 
are almost eigenvectors of the Laplacian and thus nor- 
mal modes of the Hamiltonian (uncoupled from the rest 
of the system and therefore not entangled). The progres- 
sive localization of the eigenvectors with k balances the 
growth of the correlations associated to the increase of 
k. It is the competition between these two effects what 
explains the peak for (Sk) in SF networks at moderately, 
rather that maximally, coupled nodes. 

Controlling entanglement.- The results presented 
above illustrate the particular properties that quantum 
networks introduce with respect to its classical counter- 
parts. To round off our discussion, we want to connect 
the rise-and-fall of the entanglement entropy with the 
relevant issue of the control of networks [36[ and span it 
to the quantum domain. Let us suppose that we have 
a quantum network with a rather regular topology, e.g. 
and ER graph. Our aim is to obtain the maximum en- 
tanglement with the network by adding an external node 
(our controller) as shown in Fig. 01 Two relevant ques- 
tions arise: (i) how many connections should we launch 
to the system under control? and (ii) what are the target 
nodes to be connected with the controller? 

In a random network, such as an ER graph, almost all 
the nodes are equivalent, thus we can assume that there 
are not preferred nodes to attach with, so that the num- 
ber of optimal connections remains as the most relevant 
issue to address. One would be tempted to state that the 



more connections we add the more entanglement with 
the system the controller will obtain. However, in Fig. @] 
we show that the entanglement obtained by the external 
controller increases up to a maximum at k ~ 15, while a 
further increase of its connectivity leads to the decrease 
of its entanglement entropy. The experiment illustrates 
that the amount of information we can extract from a 
quantum network via an external controller is bounded 
and its maximum is reached by means of a moderate 
number of connections with the target system. 

Conclusions.- The entropy of entanglement of nodes in 
networks of quantum oscillators reveals a non-trivial be- 
havior. The decay of the entanglement for large connec- 
tivity nodes is seen as the fingerprint of the localization 
of some of the Laplacian eigenvectors around hubs which 
turns them into normal modes of the system. This effect 
balances the increase of the entanglement with the con- 
nectivity, analogously to an area law in regular lattices, 
thus causing the rise-and-fall of the entanglement entropy 
across connectivity classes. As a byproduct, the control 
of quantum networks by maximizing the entanglement 
with an external controller is limited and does not rely 
in maximizing the number of physical connections with 
the network under study. 

In addition to their interest for the emerging field of 
quantum information networks, our results show an inter- 
esting connection with fundamental concepts of quantum 
gravity in complex spatial connectivities. In fact, the 
setup used here can be seen as the discretized version 
of real massive Klein- Gordon fields far from the usual 
Minkowsky or curved space-time situations, suggested 
from emergent gravity concepts as intermediate topolo- 
gies in the transition from a highly connected (high en- 
ergy) quantum geometric phase of the universe to the 
low energy, largely homogeneous, actual phase [37l |38|. 
Although the connection between complex space-time 
topologies and the field of network science has recently 
attracted attention [39| , the above scenario paves the way 
for the study of quantum models on complex space-time 
topologies thus fostering the synergy between these ap- 
parently unrelated fields of physics. 
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Appendix A: Entanglement entropy 

Let us start by deriving Eq. (4) in the text. Given 
g, the quantum state of the network, its associated Von 
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Neumann entropy is given by, 

S = — Tr(ghig) , 



(Al) 



where Tr accounts for the trace operation. The marginal 
entropy for the node i, Si, is obtained by replacing g by 
the reduced density matrix, 



Tr^p 



(A2) 



in (|A1[) with Ti> the partial trace. The latter means the 
trace over the complement of i (the rest of the nodes) . 

In our work, the state of the network considered is the 
ground state of the Hamiltonian [Eq. (1) in the text]: 



H network = X (P T P + X T (I + 2cL)x) , 



(A3) 



here I is the N x N identity matrix, c is the coupling 
strength between connected oscillators and L the network 
laplacian. The operators p T = (px,P2, — ,Pn) and x T = 
(xi, , X2, a; at) are the momenta and positions of nodes 
respectively, satisfying the usual commutation relations: 
[x, p T ] = ih I. We are interested in analyzing the ground 
state of the system, which is a pure state, thus having 
5 = 0. On top of that, this state is Gaussian (since the 
Hamiltonian is quadratic) so that the reduced density 
matrices ft and Si can be computed by means of the 
covariance matrix: 



-,{xiPi + piXi) 



2 {%iPi Pi%i 
„2\ 



(A4) 



where the averages are calculated via the reduced density 
matrix, g i7 as: (xf) = Tr(xfgi). 

It was Agarwal |30j who derived an explicit formula for 
the marginal entropies [Eq. (3) in the text] : 



Si — Si< 



I) log( W + I) ~ (m ~ \) log(Mi - 1) , 

(A5) 



with m = v^fKpf)- 

We are able to can find these quadratures by working 
with normal modes, i.e., those diagonalizing the potential 
energy matrix V = I + 2cL: 



SQ so that S T VS = V d , 



(A6) 



whose quadratures are those of a set of uncoupled oscil- 
lators at their individual ground state: 



{Qi, 


= {Pi) 


= 




(A7) 




(QiQj) 


= 5 *3 
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m~ 


(A8) 










(A9) 






= S ij 






f Q Pi) 


= o, 




(A10) 



with Qi = y/l + 2c\i the eigenfrequencies, and Aj the 
eigenvalues of the laplacian matrix L. The latter relation 



is obtained from simple inspection of the eigenvalue equa- 
tion Vu = il 2 u = (I + 2cL)u, so that Lu 



2c 



u = Xu. 



That is, u is an eigenvector of both V and L, with 
Q = VI + 2cA. Then we obtain: 

(x i ) = y £s ij (Q J )=0 (All) 

3 

(p i > = E^^) = ° ( A12 ) 

(jHXj) = SikSji(P k Qi) = (A13) 



(A15) 



3 3 

(fl) = E(^)>/)=E(^ 

3 3 

Finally, we arrive to the expression for the quadratures 

A*i = i x i)(Pi) - \ {xtPi+PiXi) 



«?><P?> = zE(S«) 9 («tf') a ^. (A16) 



as stated in the main text. 



Appendix B: Entanglement mean field 
approximation 

In this section we sketch the solution for the mean field 
approximation presented in the main text. The mean- 
field Hamiltonian for a node surrounded by k neighbors 
[Cf. Fig. l.a in the article]can be rewritten in matrix 
form like as: 



H = -(p*Ip + x*t>x 



whith the (k+l)-tuples: 



(Bl) 



(B2) 



note that we have named the 0-node the central one. The 
then potential reads: 





(po\ 




( X °\ 




Pi 




Xi 


p = 




, x = 






\p k ) 




\x k ) 



V 



/l + 2cfc -2c -2c ••• 
-2c 1 + 2ck 
-2c 1 + 2ck ••• 



-2c 





\ -2c 



1 + 2ckJ 



(B3) 

The equilibrium properties of (|Bip . in particular the 
Von Neumann entropy, is characterized by the eigenval- 
ues and eigenvectors of V, as explained in the previous 
section. It turns out that the spectrum of V given by 
(IB3|) can be analytically computed: 
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i) The (fe+l)x(fc+l) matrix fin (|B3j} has (ft+l)-2 
eigenvectors of the form: 

/ \ 



|Aj) 



with degenerated eigenvalues 



A,- = 1 + 2cn 



(B4) 



(B5) 



in maint text. Thus, the quadratures can be written as: 

.2 ! 



z~ 



} ^/A7 



<p3> = £ 



^ 2 + fc 2 ' 

I— ± 1 



(BH) 
(B12) 



from which the entropy is obtained. 

We finally note that z + — > —k for large enough k. 
Therefore in this limit the corresponding eigenvector ap- 
proaches to (1,0, ...,0) with frequency u> = U A + — > 
%/T-j- 2ck. Therefore, in the limit of large connectiv- 
ity the node is a normal mode and its corresponding 
marginal entropy approaches to zero. 

Appendix C: Extimation of fitting parameter k 
values 



as can be easily checked. 



ii) The other two eigenvectors are of the form, 



1 



y/z 2 + k 



1 



V 1 / 



(B6) 



where the eigenvalues and eigenvectors are found 
from the equations, 



(1 + 2ck 2 ) - kc 2 = Xz , 
-c 2 z + 1 + 2ck = A , 



(B7) 
(B8) 



with eigenvalues 



A± = - 1 2c k + yjn 2 - 2nk + k(k + 4) + k +2| 



and 



z± 



(B9) 



l r 



K - k =F \Jk 2 - 2nk + k(k + 4) (BIO) 



The latter are the only ones entering in the formula for 
the marginal entropy of the node 0, see (|A5|) and Eq. (3) 



Finally, we focus on the values of the fitting param- 
eter k used to compute the mean-field solution shown 
in Figs. 2 and 3 of the main text. The values of K for 
the different topologies under study are shown in Tab. HI 
As observed, in the case of Erdos Renyi networks with 
average degree (fc) = 2 results are missing. This is due 
to the fact that for average degrees lesser than 2.5 the 
resulting ER networks are not made by a unique con- 
nected component. As already stated in the main text, 
we observe that k only depends on the average degree 
and on the considered underlying topology. Also, as a 
further validation of our mean field hypotesis, with the 
only exception of Scale-Free networks made through the 
configurational model with (k) = 4, we notice that k 
always assumes values slightly below the average degree. 
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TABLE I. Fitting parameter k for all the different network 
topologies with respect to each pair of parameters average 
degree (k) and coupling c. 
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